On a Class of Configuration and 
Coincidence Problems 

By Z. A. MELZAK* 

(Manuscript received August 28, 1967) 

Let A and B be sets in E m where B is convex and symmetric about o. 
Let n points be taken in A and let B { be the translate of B centered at the 
i th one. Let Y be the subset oj the Cartesian product A", corresponding to 
the configurations {B x , • • • , B„) such that no more than p — 1 sets B { 
intersect, or corresponding to any similar configuration condition, expres- 
sible in purely Boolean terms. The problem of evaluating various integrals 
over Y generalizes a number of questions in queuing, telephone traffic, 
statistical mechanics of hard spheres, and so on. This article gives a complete 
solution for certain special cases, and discusses numerical {Monte Carlo) 
techniques. 

I. INTRODUCTION 

We consider here a number of problems of the following general type. 
Let A and B be two sets in the w-dimensional Euclidean space E m (m^ 1). 
B is assumed to have a center of symmetry and for any point x B(x) 
denotes the translate of B centered at x. An integer n(n ^ 2) is fixed 
and the rc-fold Cartesian product A X A X • • • X A is denoted by P. 
If u e P then u — (xi, ■•• , x„) where x { e A for t = 1, • • • , n; we shall 
be interested in the sets B(x t ), ■ ■ • , B(x„). By a configuration condition 
we shall understand a statement referring to the relative positions of 
the sets Bfa), • ■ ■ , B(x„) and describing their intersection properties 
in purely Boolean terms. 

Examples of admissible configuration conditions are: (i) the n sets 
are pairwise disjoint, (it) their intersection is empty, (Hi) their union is 
connected. A configuration condition which generalizes (i) and (n) is: an 
integer p is given (2 ^ p ^ n) and no p of the n sets intersect. Any admis- 
sible configuration condition C induces a partition of P into two disjoint 
and complementary sets Y = Y(C) and N = N(C) ; if u . = (x t , • • • , x n ) t 
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P then u t Y if and only if the condition C holds for B(z,), • • • , B(x n ). 
Finally, a function F = F(x t , • • • , x n ) is denned over P and dV denotes 
the volume element dx, • • • dx n . Our problem is to evaluate the integral 



-L 



FdV. 



In all cases to be considered the sets A, B, and Y, as well as the 
function F, will be sufficiently regular so that the questions of meas- 
urability and integrability will not arise. In fact, in most cases of 
interest B turns out to be a ball, a cube, or an w-dimensional regular 
octahedron. All these are Minkowski balls for a suitable norm />(£) = 
p ($i, . • • , $n)- We get the Euclidean ball with 

the cube with P ($) = max (&,...,£»), and the octahedron with 
i 

p(*)= El&l. 

i 

It will be therefore assumed throughout that B is a Minkowski ball. 
This amounts simply to assuming that B is a convex symmetric body. 
The precise shape of A is of no particular importance, only its con- 
tent and sufficient regularity are. 

The integrand F will be usually of some highly symmetric type 
such as 

F=l, F = flKxd, F= II K\x*-*i\)> 

t-l lS«'<j'Sn 

where / is a suitable sufficiently regular function. 

In this part of the paper we are concerned with certain special 
configuration conditions which lead to an explicit expression for J in 
terms of the so-called cluster-integrals. Later we consider a related 
expansion of the form 

j = E Ji\ j (i) 

where the parameter A measures the ratio of sizes of B to A. We shall 
take up the questions of the existence of the expansion (1) and the 
regularity of J as a function of A. 
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II. EXAMPLES 

Example 1 

Let m = 1, A is the interval [0, L] and B is the interval [0, a], n is 
any integer such that (n — \)a ^ L, the configuration condition is 
that the sets B(x l ), ■ • • , B(x n ) are disjoint, and F = 1/L. J is now the 
probability that with n points at random on the interval [0, L] no two 
points are closer than a. 

Example 2 

Let m, A, and B be as above, 

*(*i , • • • , *») = n /(*<) 

where j(x) is a probability density on A. The configuration condition 
is: p is an integer (2 ^ p ^ n) and some p-tuple of the sets B(x t ), ■ • • , 
jB(rr n ) is to have a nonempty intersection. Here we have the following 
interpretation: [0, L] is a basic time interval and n events occur during 
that time. Each event occurs independently of the others with the 
probability density j(x). A p-fold coincidence is defined to be the com- 
pound event arising when some p events occur closely together — on a 
time-interval of length a. Now J is the probability that a p-fold coinci- 
dence occurs. 

The above examples show that problems of our type might be of 
interest in queuing theory, telephone traffic, the theory of particle 
counters, and in similar areas. The next example is a scattering 
problem for a random linear array of n identical isotropic point- 
scatters, no two of which can be too close together. 

Example 8 

Let m, A, B, and C be as in example 1. We suppose that the wave- 
length is 2ir and that L is an integral multiple of it. Aside from propor- 
tionality factors the signal scattered by the array is the vector (|, 77) 
where 

£ = £_, cos x { , 7} = ^ sm x < • 

1 1 

We are here interested in the probability P(u, v) that 

tt ^ I :£ u + du and v ^ 77 ^ v + dv. 
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The Markov method 1 gives 

P(u,v) = (27r)-V.r f f X e- i(T " + " ) B(r,s)drds 

J— oo J — oo 

where 

B(r,s) = [ exp i[r £ cos *'. + s £ sins, J dV. 
J Ya - \ i ' /J 

J.4 and 7.4 are the integral and the region of example 1, respectively. 

Therefore the spectrum B(r, s) is obtained in the form of our integral 

J if we take 

F= fife). /(*) -•«'- — — , . 

i 

When a = then P{u, v) reduces to the probability density for the 
isotropic plane random walk of n unit displacements in arbitrary 
directions. 

Example 4 

Let m — 3, let 4 be any large and sufficiently regular portion of 
space, and let B be the ball of radius a. The configuration condition 
is that no two sets B (x t ) and B (x } ) overlap. There is a suitable given 
function f (x) and 

F(x lt ■■■ ,x n ) = n •- <, " i -" ) . 

lSi<i&n 

Now, aside from some simple normalization factors, J is the so-called 
partition function for a hard-sphere model of idealized gas with inter- 
molecular potential <p and the hard core radius a. 2 

The knowledge of J is here of considerable importance in statistical 
mechanics and a great deal of work has been done on the subject of 
evaluating J in the form (1) which is closely associated with the so- 
called virial expansion. 

III. A SPECIAL CASE 

The method to be used involves certain dissections of Cartesian 
products together with the inclusion-exclusion principle of combina- 
torics. 3 As an illustration and an introduction to the more complex 
examples which follow, we consider here at some length example 1 
of the previous section. The material is taken from Ref. 4, where some 
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further details can be found. The well-known solution 5 - ° is here 

./ = J(n, a, L) = [1 - (n - l)a/L] n (2) 

and it may be obtained analytically as follows. 

Let the coordinates of the n points be x 1 , ■ ■ • , x n ; these can be ordered 
in n ! ways. Suppose that ^ x t ^ x 2 ^ • • • ^ x„ ^ L ; the conditions 
of the problem are satisfied if and only if 

^ x t ^ Xi - a ^ x 3 - 2a 

^ • • • £ Xn - (/i - l)o < L - (n - l)a. (3) 

Let iji = x t — (i'-l)a (i = 1, . . . , n), then the probability that (3) 
holds is L _n times the volume of the region in E" consisting of the 
points y = {y lt . . . , y n ) for which 

^ y x ^ ?y 2 ^ • • • ^ !/„ ^ L - (n - l)o. 

The volume in question is [L— (?? — l)a n ]/n\; since there are ?i! equi- 
probable orderings we get (2) at once. 

Consider next an alternative geometrical proof of (2), which is 
considerably more complicated, but leads to useful generalizations 
and gives some additional insight. 

First, let n = 2. The sample space of pairs (x x , x 2 )(0 ^ x t , x 2 ^ L) 
is the square Q of side-length L, lying in the first quadrant of E 2 and 
containing the origin as a vertex. Let D be the diagonal of Q through 
the origin and draw the two lines parallel to D at the distance 2~*a 
from it. The hexagonal subset of Q contained between those two lines 
is the sample space of the forbidden configurations with |x, — x 2 \ ^ a. 
The remainder of the square Q consists of two congruent triangles which 
can be moved together so as to form a square Qi , of side-length L — a. 
By the randomness assumption J (2, a, L) is the ratio of the areas of 
Qi and Q which yields (2) for n = 2. 

The case of arbitrary n is handled similarly. In E" we take a 
Cartesian coordinate system with the n axes X lf . . . , X„. The n- 
dimensional cube 

H - \(x lt • • • , .t„): ^ .r,. g L, % = 1, ••• , n) 

is then the sample space of all /^-tuples of points on the segment [0, L]. 
Let 7, be the interval [0, L] on the X { - axis. In the two-dimensional 
square face Q ti = 7, X 7, of 77 let D,, be the diagonal through the 
origin and let 77,, be the hexagonal subset of Q ti consisting of all points 
no further from D it than 2~*a. Let £,,- be the Cartesian product of 77,, 
with all the I k 's for which k ?± i and k ?± j. 
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Sa is now the sample space of the configurations which are forbidden 
on account of too close approach of the points x { and x,- for the chosen 
indices i and ;: \x { — rc,| ^ o. The sample space Y of the allowed con- 
figurations is therefore the set 

H - U B» ■ 

lSt<lSn 

When the (z\ paradiagonal slabs £,-,■ , based on the paradiagonal sets 

En , are removed from H, the remainder of the cube H consists of n! 
congruent simplexes which can be reassembled by suitable translations 
so as to form a smaller cube H x of sidelength L — (n — l)a. By the 
randomness assumption Jin, a, L) is the ratio of the volumes of the 
cubes H ! and H, and so (2) is proved again. 

The above procedure works on account of a lucky geometrical accident 
of the fitting of n! simplexes. If A and B were some other, ra-dimensional, 
sets, we could still form the paradiagonal sets and slabs and we could 
attempt to find the volume of the union U S ti of all the paradiagonal 
slabs. This is essentially what is done in the next section by means of the 
inclusion-exclusion principle . 

IV. SIMPLE COINCIDENCE WITH SEPARABLE INTEGRAND 

In this section we are concerned with a configuration condition cor- 
responding to simple coincidence : u = (x x , • ■ • , X n ) e Y if and only if 
for some i and j B{x { ) and B{x { ) intersect. Subject to general restric- 
tions, A,B,m, and n are arbitrary. We let N = u] and we form the N 
paradiagonal sets 

and the N paradiagonal slabs 

S„ = [(x t , ■■■ ,xJ:BWf\BM *4>\. 

Let the slabs be enumerated by a single index as {S k }, k = 1, . . . , N. 
Then an application of the inclusion-exclusion principle gives 

/= [FdV= ZC-irf E-Z f FdV\ (4) 

Jy r -l Llfi*.<*,<---<A:rS-V J S k ,r\- • -ASl, J 

= 2(-irx. 
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With the general integrand no further elaboration of (4) is possible. 
Suppose now that F has the separable form 

j-n/<*). (5) 

With the double-index enumeration of the iSy's the first term K x can 
be written as 

k* = z f n aw dv 

and since all the N paradiagonal sets are congruent, we have 

Ki = N iL /(x) dx ) 2 / K x M x ^ ^ dx * • 

For reasons which will be clear shortly we write 
N u = N, I /(as) dx = J , f f(x l )f(x 2 ) dx, dx 2 = /„ (6a) 



so that 



Kx = NuJT'Jn • (6b) 

Similarly, the second term K 2 in (4) is 

(•t.i'i) («t.j.> • l Si 1 i 1 r,Si,,- t i 

where the summation extends over all distinct pairs (i, , ; x ), (» a , ; 2 ) 
such that 1 ^ t| < J, ^ n, 1 ^ i 2 < j 2 ^ n; no regard is paid to the 
order of pairs; [(1, 2), (3, 4) is the same as (3, 4), (1, 2)] so that there 
are exactly 

(?) 

such pairs of pairs. There are two types of these: N 21 pairs like (1, 2), 
(3, 4) with all four indices different, and N 22 pairs like (1, 2), (1, 3) with 
one shared index. By a simple calculation 

N 2l = n(n - \){n - 2)(n - 3)/8, N 22 = n(n - l)(n - 2)/2, 

(7a) 



N n + N 22 



-(?) 
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and in analogy to (6) we set 

Jn = / PI 1( x <) dx i dx * dx * dx * 

= \f /(.rO/(x 2 ) dx x dx 2 j = Jl (7b) 

j 22 = I n K x i) dx > dx * dx * 

■';/,, mil, l 

so that 

K 2 = N 2 iJT*Jli + N 22 JT*J 22 . (7c) 

The main purpose of this section is to develop formulae, analogous 
to (6) and (7), for the general term K r of (4). The principal dif- 
ficulty here is that in passing from the single-index formula for K r 

K r = s • • • E ( n Kxi) dv m 

to the double-index formula 

K r = E ••• E ( ILK**)dV (8b) 

(...i.) (ir.fr) •'S( t y i rV"r\S,- r ; I . 1-1 

we need an adequate description of the different types of r-tuples 
of pairs of indices occurring in (4) , together with a hold on the range 
of summation in (8b). For instance, with r = 2 there are two such 
types, illustrated by (1, 2), (3, 4) and (1, 2), (1, 3). With r = 3 there 
are five types of index-sharing in triples of pairs: 

(1, 2), (3, 4), (5, 6); (1, 2), (1, 3), (4, 5); (1, 2), (2, 3), (3, 4); 

(1, 2), (1,3), (1,4) ; (1,2), (1,3) (2, 3); (9) 

We may therefore expect that the formula for r = 3, analogous 
to (7c) for r = 2, will have five terms rather than two. The number 
of such types grows very rapidly with r, and as an aid we introduce 
certain graphs associated with the terms of (8). These graphs reflect 
completely the intersection properties of the sets B{x x ), . . . , B(x n ). 
For r = 3 there are five such graphs corresponding to the five types 
enumerated in (9). These are given in Fig. 1 together with the cor- 
responding .^-configurations. (It is, of course, assumed that n > 2.) 

Each graph is of the following kind: 

(i) No vertex is isolated. 
(ii) No pair of vertices is connected by more than one edge. 
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Fig. 1 — Coincidence graphs, r = 3. 

(Hi) No edge connects a vertex to itself. 
(iv) There are exactly r edges. 
(v) There are exactly v vertices. 

One further, and crucial, condition is added: 

(vi) If the v vertices are enumerated in some order then there exists 
a configuration of v translates B t , • • ■ , B, of B, such that B { 
and B, intersect if and only if the I th and the ; th vertices are 
connected by an edge. 

For the sake of convenience we make here the following conven- 
tion: two convex m-dimensional bodies will be said to intersect only 
if their intersection is itself w-dimensional, otherwise they are to be 
regarded as disjoint. The reason for this is that we are interested in 
purely metric properties: the intersections of such sets serve as do- 
mains of integration for well-behaved functions in E' n . 

A graph satisfying conditions i through vi will be called a [B, r, v)- 
graph, one satisfying i through iv and vi a (B, r) -graph, and one 
satisfying i through Hi and vi a B-graph. It must be emphasized that 
the condition vi is not of the usual graph-theoretic kind and it pre- 
vents many graphs from being B-graphs. For instance, let m = 2 
and let B be a circular disk. Since a disk in E- cannot intersect six 
congruent pairwise disjoint disks, the graphs of Fig. 2 are not B- 
graphs. 

The proof of the above assertion for the graph of Fig. 2b is obtained 
by showing that here the "extreme" configuration is that of Fig. 3. 

Similarly, when m = 2 and B is a square then B cannot intersect 
five pairwise disjoint translates of itself (for each translate contains 
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05 




Fig. 2 — Graphs which are not B-graphs. (B is a disk.) 



a vertex of B) so that the graph of Fig. 4a is not a B-graph. On the 
other hand, the graph of Fig. 4b, which corresponds to that of Fig. 
2b, is a B-graph as shown by the configuration of Fig. 4c. 

Returning to the evaluation of K r , we start with (8b). Summation 
there extends over all the 

G) 

2 

distinct r-tuples of pairs of indices where for each pair (i, , j.)l =a i, < 
j, ^ n; r-tuples differing only in the order of pairs are not considered 
distinct. We can now associate the terms of (8b) in a 1 : 1 fashion with 
the distinct (B, r)-graphs on some n vertices w x , • • • , w n . Given a 
B-graph G let 

S(G) = n S<* . dO) 

where the intersection is taken over all pairs {i, j) for which w t is con- 
nected to Wj by an edge in G. Then (8b) may be written as 



K r = z [ n m dv, 



(11) 




Fig. 3 — An extreme B-configuration. 
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Fig. 4 — Configurations when B is a square. 
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the summation running over all distinct (B, r) -graphs on n vertices. 
Let v (G) denote the number of vertices of G and C(G) a connected 
component of G. Since the integrand in (11) is completely separable, 
the integral over S(G) splits into a product of integrals over the con- 
nected components and we get 



K r = Z Jo-' (G) II J[C(G)). 



(12) 



Here J[C(G)] is an integral over the connected component and the 
product is taken over all such components of G. Two examples of 
integrals J[C(G)] are given in (7b). Owing to the congruence of all 
the paradiagonal slabs and the form of the integrand, it is not neces- 
sary to sum in (12), over all (B, r) -graphs on the vertices Wi, . . . , 
w„, but only over their types. 

Suppose that there are exactly t = t(r) types of such graphs and 
let Gj be any one of the ;' th type; let also N rj (n) be the number of 
different (B, r) -graphs on the vertices u>i, . . . , w n , of the j th type. 
Then (12) becomes 



Mr) 



K T = ZN ri (n)r -" ao II J[C(G,)]. 



(13) 



Thus the problem of evaluating J has been reduced through (4) 
and (13) to: the geometrical problem of determining the types of 
(B, r) -graphs, the combinatorial problem of calculating the coef- 
ficients N r} (n), and the analytical problem of evaluating the cluster- 
integrals over the connected (B, r) -graphs. 

V. MULTIPLE COINCIDENCE WITH SEPARABLE INTEGRAND 

Formulae analogous to those of the previous section will now be ob- 
tained for the case of p-tuple coincidence. Subject to general conditions, 
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A, B,n, and m are arbitrary and F is of the separable form (5). An in- 
teger p is fixed (2 ^ p ^ n) and the configuration condition is: u = 
(xi , • • ■ , x n ) e Y if and only if there are p indices i x , ■ ■ • , i v (l ^ *i < 
< *'p ^ n) such that 

j = l 

We observe here our convention that the intersection must be it- 
self Hi-dimensional. We introduce the analogs of paradiagonal sets 
and slabs: 

//,,...,„ = |(z„ , • • • , *J : fj B(x<.) * 4>J , 

we let Af = I"), and we re-enumerate the M sets £,,...,,, with a single 
index k as {S k \, 1 ^ k ^ M. Then we get a formula analogous to (4): 

/- [ Fdv = ec-We--- z f ^ F ] 

- Z(-i) r+1 ^- • < 14) 

As in (6a) we let 
M u = M, [ f(x) dx = J , I II /(»<) ^i • • • dx > B J " • 

J. -I J Hi,... P 1 

to get 

In terms of p-tuple indices the second term C/ 2 of (14) is 



The summation extends over the f J distinct pairs of p-tuples. We have 

now p types of such pairs, depending on the number of shared indices, 
which may be 0, 1, • ■ • , or p - 1. Let M 2i be the number of p-tuple 
pairs of type j (that is, with y — 1 indices shared) and put 

,, 2p-i+l 

JH>,...pr\ii P +,-j..., P -i+> i-i 
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then 

Observe that the integral J 21 splits into a product: J 2i = J 2 U . To 
get an expression for arbitrary U r we introduce a higher-dimensional 
equivalent of E-graphs. Let X be a regular simplex in E n ~ x on the vertices 
w r , • • ■ , w n . On account of properties i through Hi listed in Section IV, 
a (B> r)-graph is simply a set of certain r edges (or one-dimensional 
faces) of X. A d-dimensional hypergraph G will be just a set of some of 

the ( J d-dimensional faces of X. This takes care of properties 

i through Hi. When there are r such faces in G we shall speak of an 
(r)-hypergraph and when these faces comprise between them v vertices 
of X, G will be called an (r, p)-hypergraph. 

An equivalent of the important condition (vi) is very naturally 
obtained : there is a /^-configuration of v translates B x , • • • , B, of B, 
such that any d + 1 of them, say, B it ,..., fi+t intersect if and only if 
w <i > " * ' , w ii+l are the vertices of a d-dimensional face of X included 
in G. Components, types, and so on, for (B, r)-hypergraphs are defined 
in the same way as before. For instance, a hypergraph G is connected 
if no plane disjoint from it can strictly separate some of its d-faces from 
others. All quantities such as C(G) and v{G) have the same meaning 
as before. Let t = t (r, d) be the number of different types of (B, r)- 
hypergraphs, let (?, be any one hypergraph of the j th type, and let 
M d rj {n) be the number of different (B, r)-hypergraphs of the j th type on 
the n vertices. Then, proceeding as before, we get the equivalent of (13) : 

u. = "zfVrW"-'' '" n J[c(G t )]. (15) 

CiGj) 

VI. SOME COMBINATORIAL PROPERTIES OF B-GRAPHS AND B-HYPERGRAPHS 

Let ip(r) and ^(r) be the smallest and the largest number of ver- 
tices, respectively, in a (B, r, v) -graph G. From conditions i through 
Hi we have at once xf/(r) = 2r. G is then minimally connected with r 
components (Fig. 5a). Suppose that r is a triangular number: r = 
s(s-l)/2; there is then a complete graph on s vertices which is 
clearly a B-graph for any B, so that s = v. If r is not a triangular 
number let t(t-l) < 2r < t(t+\) and put e = r-t(t-l)/2. 

Let G be the complete graph on t vertices. For the corresponding 
.B-connguration we may assume that the translates B t , • • • , B t of B 
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r-20— o r-20— — o r-20— o 

(a) (b) (c) 

Fig. 5 — (B, r, v) — graphs with high v. 

have an interior point in common. We check that e < t and that B lt - ■ • , 
B, may be arranged so that a point z s f\'Bi can be strictly separated 
from UJ + i Bi ; by a plane P. Let B, +1 be a translate of £ which contains 
z and lies strictly on the same side of P as z. Then the resulting B- 
configuration B x , • • • , B t+1 corresponds to a (B, r, v)-graph G with 
„ = 1 _j_ i_ This (? may be said to be a maximally connected (B, r, v)- 
graph. We have now 

vKr) - 2r, <p(r) = min {j : j f fc [1 + (1 + 8r)*]/2}. 

Similarly, let ^(r, d) and ^(r, d) be the corresponding minimum and 
maximum of v for a (£, r, y) -hypergraph. Then clearly f(r, d) = (d + l)r. 

To determine <p(r, d) we suppose first that r = ( , , j. There is then 

a complete hypergraph on s vertices, consisting of all the d-dimensional 
faces of an (s - l)-dimensional simplex. This hypergraph is a 5-hyper- 
graph for any B and so v = s. If 

Ux)<'<o:i) 

we proceed as before and find that v = t + 1. Hence 
yp{r, d) = (d + l)r, <p(r, d) = min [j : 

j ^ largest pos. root of x (x - 1) • • • (x - d) = (d + l)!r|. 
The bounds p(r) and i^(r) lead us to the possibility of a combina- 
torial identity 
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and its relation to J3-graphs and the numbers N r j{n). For instance, 
we find for r = 2 

?)■•© + -©■ 

2 J is the total number N 2l + iV 2 2 of graphs of (7a) and 

3© = *,,, 3(™) = ^. 

To prove the validity of an expansion like (16) for all n observe 
that the left-hand side is a polynomial in n of degree 2r = \p(r) so 
that 

Further, A rk = for k < <p(r), for we substitute successively n = 0, 
1, • • • , <p(r) — 1 in (16) and recall that ( p j = for p < r. By expanding 
both sides of (16) in powers of n and comparing the coefficients we find 
A r%r = (2r)!/2'r!, A r2r _ t = (2r - l)\/7T l (r - 2)!, 
A r2r _ 2 = (2r - 2)!(3r - l)/8.2 r " 1 (r - 3)! 
and so on. Therefore (16) may be written as 

(?) - (»W(2V!) + (n) 2r ^/[2 r -\r - 2)!] 

+ (?i) 2r - 2 /[3.2 r -\r - 3)! (3r - 1)] + ■ • • (17) 

2r-#>(r) 
i-O 

(n) p stands for »(n— ) . . . (n—p+1). 

The denominators D t have the following interpretation. Consider 
first the {B, r, 2r) -graph of Fig. 5a. The 2r vertices can be chosen out 
of Wi, . . . , w n in (n) 2r ways. We define the symmetry number for 
a (B, r, v) -graph to be the number of ways in which its vertices can 
be labelled with integers 1, 2, . . . , V, all of which ways are to corre- 
spond to the same ^-configuration. Here the symmetry number is 2 r r!, 
as there are 2 r ways of permuting the labels on the two vertices of a 
component and r! ways of permuting their components. This leads 
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us to the first term (») 2r /2 r r! in (17) which is precisely the number 
AT rl (n) of (13) provided that we consider G, in (13) to be of the type 
of Fig. 5a. 

Similarly, for the (B,r,2r-1) -graph of Fig. 5b we find the sym- 
metry number to be 2 r ~ 1 (r— 2)!. The number of ways to choose the 
2r— 1 vertices is (n) 2r -i and so we get the second term (n) 2 r-i/2 r_1 
(r— 2) ! of (17). The situation gets somewhat more complicated for 
the (B,r,2r-2) -graphs. Here we have three types instead of one, 
illustrated in Fig. 5c. The 2r-2 vertices can be selected in (n) 2r _ 2 
ways, the symmetry numbers for the three types are 

2 r ~\r - 3)!, 3.2 r_2 (r - 3)!, and 2 r_1 (r - 4)!. (18) 

Therefore, the corresponding numbers of graphs, say N r ^(n), N rA {n), 
N r5 {n) are 

(r*) 2r _ 2 /[2 r - 2 (r - 3)!], (n) 2r _ 2 /[3.2- 2 (r - 3)!], (tt) 2r _ 2 /[2 r-1 (r - 4)!] 

and their sum is precisely the third term of (17). The corresponding 
denominator D 2 is therefore three times the harmonic mean of the 
three symmetry numbers in (18). 

Thus the first few terms of (17) give the total numbers 

i 

of (£,r,i>) -graphs for v = 2r,2r—l, and so on. However, this pleas- 
ing circumstance breaks down as soon as we reach the smallest term 
i for which one of the types of graphs in question is not a B-graph. 
For the case m = 2, B a circular disk, this occurs for i = 7 and the 
graph in question is then that of Fig. 2a together with other components 
containing one edge each. When B is a square the graph of Fig. 4a shows 
that the breakdown occurs for i = 6. On the other hand, the quantity 
(n) 2 r-i/Di from (17) always provides an upper bound for the sum 
2 Nn(n), the summation extending over all types ; of (B, r, 2r — i)- 
graphs. 
The explicit form of (16) is 

where q = y(r) and 
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We prove (2) by induction on k. For k = q, (20) holds, suppose it to be 
proved for k ^ q + — 1. Let n = q + s in (19), then 

which by the induction hypothesis may be written as 

^,..=(™)-Ei(-.)'(':f ri) 

\ r I frT. TT n \q + i/\ r 

In the double sum we may sum first over those terms for which the 
difference u = i — j is constant, then over u. In this way one gets 

which after some simple algebra becomes (20) with k — q+s. This 
completes the induction and the proof of (20). 

Some combinatorial identities may be obtained from the above. 
For example, we know that A r2r = (2r) !/2'r!. Hence, on putting k = 2r 
in (20), we get 

Z(-i)'0( r ;' ) ) = (2 '' )!/2V! - (21) 

Similarly, with k = 2r — 1 and k = 2r — 2 we get 

2r- <i-l 



"l (-l)'f r • l )(f f 'p) = ^ ~ 1)!/2r_I(r " 2)! (22) 



and 



2r-a-2 



= [(2r - 2)! (r - l/3)]/2 r -(r - 3)!. (23) 

For hypergraphs we have the identity 

where q = f (r,d). The explicit expression for the coefficients A rk (d) 
can be found in the same way as (20) : 
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Some of the higher coefficients A r dn A r dr _i , ... can be evaluated by 
comparing the powers of n in (24) : 

A rdT (d) = (dr)\/r\{dW, A rdr _M = (dr)\d(r - l)/2.r\(d\) r 
and so on, so that by putting k = rd and k = rd - 1 in (25) we get 

23 (-l)'( (dV) ) = (dr)\/r\m r (26) 

,--0 \ 7 / 



and 

dr-Q-l 



ff (_l)'( rfr " ^(C'-r 1 )^ = (dr)]d(r _ i)/2. r !(d!)'. 
,=o \ j / \ i l 



(27) 



The coefficients A kr (d) have the same interpretation with hyper- 
graphs as the A kr have with ordinary graphs, and they refer to 
symmetry numbers. 

VII. SIMPLE COINCIDENCE IN A CUBE 

We consider here the problem of evaluating the probability P(n,a,L) 
that when n points are taken at random (uniform distribution) in a 
three-dimensional cube of edge-length L, then no two points are closer 
than a. The problem occurs in deriving the van der Waals equation 
from a primitive hard-sphere gas model. See, for instance, Ref. 2, 
where the problem is termed "very difficult" and the crude (though 
sufficient) approximation 

P(n, a, L) ^ n (1 - 47rja 3 /3L 3 ) ^ 1 - 2 7 m 2 (a/L) 3 /3 (28) 

is used. 

From our formulation we find that 

L 3n [l - P(n, a, L)] 

is the J integral for the case m = 3, A is a cube of volume L 3 , B a ball 
of radius a/2, and the configuration condition is that not all sets B (Xi) 
be disjoint; in other words, a simple coincidence. Therefore by (4), 
(13), and an inspection of Fig. 1 we have 

L 3n [l - P(n, a, L)) = N n L 3 "- 6 I u - (N 2l L 3 "- 12 I 21 + N 22 L 3 "- 9 I 22 ) 

+ (Ar 31 L 3n - 18 J 31 + JV 32 L 3n - ,5 J 32 + N 33 L 3n ~ 12 h 3 

+ AT 34 L 3n - ,2 / 34 + N 35 L 3n -% B ) (29) 
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where the integrals In, hi, ■ . . can be symbolically represented as 
follows 

mi — / , I21 ~ I — I"n , la* — I t I31 = I = /?i , 

~ ~ v i (30) 

I32 — / = /n/22 . 

A~ 

To obtain an explicit Cartesian expression for an integral, I, we 
consider its signature graph G which is a (B,r,v) -graph. If the v 
vertices are enumerated as 1,2, . . . , v in an arbitrary order then J 
becomes a 3t>-tuple integral 

I = j ... f dn-'dr, (31a) 

where r< is the vector (n\, j/ ( , 2 4 ), dr< stands for dx t dy t dz iy and the 
region of integration Rj is given by 3v + r inequalities : 

£*,££, £ /y, gI,0^^L, (,'-!,..., y ), (31b) 

|r, — r,-| 2 ^ a 2 if the i th and the ;' th vertices are connected (31c) 

in G by an edge. 

Further, such an integral occurs in (29) with the multiplier N rj L sn ~ zv 
where N rj is the number of distinct graphs on n vertices, which are of 
the same type as G. Together with each such integral I = I pq we may 
also consider the corresponding integral K pg given by 



K vq = J • ■ ■ J dr x • • • dr, , 



where the region Q, is given by the (v 2 + 5y)/2 inequalities (31b), (31c) 
and 

|r< — Tj\ 2 ^ a 2 if the i li> and the j th vertices are not connected in G 

(31d) 

by an edge. 

It turns out that the J integrals are expressible in terms of the K 
integrals, and conversely. For instance, consider the K integral with 
the signature graph which has four vertices 1, 2, 3, and 4, and edges 
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12, 13, 23, and 24. We write it in a self-explanatory terminology as 

(12)(13)(23)(24)[1 - (14)][1 - (34)] 
and multiply this out to get 
(12) (13) (23) (24) - (12)(13)(23)(24)(14) - (12) (13) (23) (24) (34) 

+ (12) (13) (23) (24) (14) (24) (34) 

which yields at once a representation of K as a sum of four /-integrals. 
The first integral I u is sixtuple and can be reduced to an iterated 
integral as follows: 

■■■ \ dXi. dx 2 dyi dy 2 dzx dz 2 (32) 

rn 2 = m 4 = m 6 = 0, M 2 = M 4 = M b = L 



where 

and 



m! = max {0, x 2 - [a — {y x - y 2 ) 2 — (z, - z 2 ) 2 ]*}, 
m l = min {L, x 2 + [a 2 - {y^ - y 2 ) 2 - {z, - z 2 ) 2 ]*}, 
m 3 = max {0, y 2 - [a 2 ){z 1 - zj 2 ]*}, 
Mi = min {L, y 2 + [a 2 - (z, - z 2 ) 2 ]*}, 
m 6 = max {0, z 2 — a\, 

M B = min {L, z 2 + a}. 

This arrangement of the limits of integration corresponds to taking 
two balls of radii a/2 and centers {x x , yi, Zi) and (x 2 , y 2 , z 2 ), and 
letting the center of the first ball move freely over the cube while the 
coordinates of the second center vary so that the balls intersect. 
Accordingly, In has a simple probabilistic interpretation: In = L 6 [l — 
P(2,a,L)], where P(2,a,L) is the probability that two points taken 
at random in the cube of edge-length L are no nearer than a. Similar 
probabilistic interpretation holds for any other K integral. If G is its 
(B,r,v) -graph then K is L 3V times the probability that when v balls 
of radius a/2 are taken with their centers at random in the cube, then 
the balls are in the configuration of G (so that two of them intersect 
if and only if the corresponding vertices of G are connected by an 
edge). 

We evaluate now the integral (32) subject to the condition a ^ L. 
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Integration with respect to x x and x 2 gives 

L 2 - [max (0, L - D)f (33) 

where 

D 2 = a 2 - ( Vl - i h ) 2 - ( Zl - z 2 ) 2 . 

Since o|Lwe have D ^ L and therefore (33) is 2LD — D 2 . Integrating 
this with respect to y t and y 2 we get first, on putting y x — y a = u, 

/ [2L(6 2 - u*)* - (b 2 - u 2 )] du dy 2 

n J-m 

where 

b 2 = a 2 — (2, — 2 2 ) 2 , wi = min (?/ 2 , b), M = min (L — y 2 , b). 
Again, a ^ L implies b ^ L and the double integral is therefore 

ttL% 2 - 8L6 3 /3 + fo 4 /2. 
Finally, integrating with respect to z x and z 2 we get 

7„ = 47ra 3 L 3 /3 - 37ra 4 Z,72 + Sa 5 L/5 - a°/Q, g a g L. (34) 

There are two more forms of I u , corresponding to the ranges L ^ 
a ^ 2*L and 2*L ^ a ^ 3*L, but they do not appear to be expressible 
in terms of elementary or standard transcendental functions. It may 
be observed that the leading term in (34) is the product of the volumes 
of the cube and the ball of radius a. 

To get a better approximation to P{n,a,L) than (28), we examine 
(29) and find that for small a every integral I ih beyond J n , is (a 8 ). 
Therefore 

P(n, a, L) 

= 1 - (£)[4*/3(d/L)> - 37r/2(a/L) 4 + 8/5(a/L) 5 ] + 0[(a/L) 6 ]. (35) 

It is possible to find the exact limit of P{n,a,L) as 

n -> oo, a -> 0, (47r/3)(n72)(o 3 /L 3 ) -♦ b. 
For we have then P(n,a,L) = P(b) and 

1 - P(b) = Nnln/L* - N 2l I 21 /L 12 + N 31 I 3l /L ia - • • • 
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and 



y M s(®)ss»W 



/*, = (In)*. 

This amounts to neglecting all graphs other than the "principal" one, 
for each k, that is, the one corresponding to the configuration of Fig. 
5a. Hence 

1 - P(b) - S (-l)' +1 [(47r/3)(n 2 /2)(a/L) 3 ]7i! = 1 - «" 6 

so that 

P(6) = e"\ (36) 

VIH. NUMERICAL EVALUATION OF THE I-INTEGRALS 

Since no I integral beyond In appears to be explicitly evaluable in 
terms of standard functions, the possibility was investigated of com- 
puting those integrals numerically by the Monte Carlo method. The 
first set of trial calculations was performed on I u itself, in order to be 
able to compare the results with the known true value. We assume 
as before that a < L and we put L = 1 (homogeneity!) to get 

I n (a) = 4.18S8a 3 -4.7129a 4 +1.6000a 8 -0.1667a 6 , £ a £ 1. 

We now choose a suitable integer M and set the value of a at 1/M. 
Next, two points pi(x lt y lf Zi) and p 2 (z 2 , 2/2, z 2 ) are taken at random 
in the unit cube by choosing each coordinate to be a random number 
from the rectangular distribution on [0, 1]. Such pairs of random 
points are selected N times; suppose that in Ni of them the distance 
between the random points does not exceed 1/M, then the quotient 
Ni/N is taken as the Monte-Carlo approximation to In (1/M). Then 
the whole procedure is repeated with 1/M replaced by 2/M, 3/M, and 
so on, until the value 3 1/2 is passed. The whole calculation will be 
referred to as an N by M Monte Carlo run. 

In the first set of trial computations N by M Monte Carlo runs 
were executed for various values of N and M, and in each case a 
least-squares fit was done on these data by a polynomial of the form 

tv. 

1-3 

The results are shown in Table 1 . 
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Table I — First Trial Computations 
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True value of A,- 


4.1888 


-4.7129 


1.6000 


-0.1667 


1st 1000 by 20 run 


3.1873 


-1.0742 


-2.5488 


1.3448 


2nd 1000 by 20 run 


3.3296 


-1.6727 


-1.9088 


1 . 1584 


10000 by 20 run 


4.4765 


-5.9918 


3.4012 


-0.9760 


1000 by 200 run 


4.3008 


-5.2689 


2.4437 


-0.5641 


10000 by 200 run 


4.1974 


-4.7337 


1.6358 


-0.1911 


100000 by 20 run 


4.1546 


-4.5615 


1.4043 


-0.0879 



It appears from this polynomial that very long and large runs are 
necessary to determine the coefficients with fair accuracy. However, 
the values of the integral itself can be computed quite well. To check 
this we have computed the standard deviations, both for the Monte 
Carlo data, from 

a\ = 1/M Z [(N,/N) - /„(i/M)] 2 
and for the least squares fit from 



where 



a\ = 1/M Z [In(j/M) - /n(i/M)] 2 



/„(«) = Z A ,11* 



is the least-squares fit to J u . The results are shown in Table 2. 

As a compromise between accuracy and length of the Monte Carlo 
run, the values N = 10000 and M = 20 were selected. In this way 
there were computed the two integrals 7 3 i and J 32 corresponding to 
the two (B,r,S) -graphs, the six integrals J 41 , . . . , 7 46 corresponding 
to the six (5,r,4) -graphs, and the 21 integrals hi, ... , h 2 i correspond- 
ing to the 21 (B,r ,5) -graphs. The first two series are shown in Figs. 
6 and 7. The programming was quite simple and no details need be 
given. The total time taken up on the CDC 6600 computer was about 
one hour; this, however, includes a lot of trial runs and tests. 

Table II — Standard Deviations 



N 


Monte 
Carlo vi 


Least 
Squares at 


1000 

10000 

100000 


0.01154 
0.00257 
. 000922 


0.00753 
0.00183 
0.000554 
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Fig. 6 — Cluster integrals for v = 3. 
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Fig. 7 — Cluster integrals for v = 4. 
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To sum up, it appears that numerical computation of ./-type inte- 
grals is quite feasible, with the help of an automatic computer, to 
fairly good approximation. One well known advantage of the Monte 
Carlo method of evaluating multiple integrals was clearly brought 
out; namely, its relative independence of the dimension. 
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